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ABSTRACT 

We present Spitzer IRS data on 19 asymptotic giant branch (AGB) stars in the Large Magellanic 
Cloud, complementing existing published data sets of carbon stars in both Magellanic Clouds and the 
Milky Way, to investigate the effects of metallicity on dust and molecular spectral features arising 
from the circumstellar envelope. We find that the C2H2 P and R branches at 7.5 micron are affected 
by dust dilution at higher mass-loss rates - albeit to a lesser extent for sources in the Magellanic 
Clouds, compared to the Milky Way - while the narrow 13.7 micron C2H2 Q branch only shows the 
effect of dust dilution at low mass-loss rates. A strong metallicity dependence is not observed for the 
Q branch. Independent of metallicity, we also provide an explanation for the observed shifts in the 
central wavelength of the SiC emission feature, as we show that these are largely caused by molecular 
band absorption on either side of the dust emission feature, dominating over shifts in the central 
wavelength caused by self-absorption. We have devised a method to study the dust condensation 
history in carbon-rich AGB stars in different metallicity environments, by measuring the strength of 
the 11.3 [Mm SiC and 30 fiin MgS features with respect to the continuum, as a function of mass- loss 
rate. With this method, it is possible to distinguish in what order SiC and graphite condense, which 
is believed to be sensitive to the metallicity, prior to the eventual deposit of the MgS layer. 
Subject headings: Magellanic Clouds - stars: AGB - stars: carbon - stars: mass loss - stars: infrared 
dust - circumstellar matter 



1. INTRODUCTION 

Stars with a main-sequence mass of ~l-8 M undergo 
an asymptotic giant branch (AGB) phase which plays an 
important role in the recycling of stellar material into 
the ISM. D uring the AGB ph ase, a significant fraction 
(up to 90%: IWeidemannlll987f ) of the stellar mass is lost 
in a gradual stellar wind. Initially, this wind results in a 
mass-loss rate of 10 -7 M yr _1 and eventually evolves 
towards a sup erwind phase with M ~ 10~ 4 M yr _1 
(Habing 1996). Given the shape of the initial mass 
function, AGB stars contribute a substantial amount of 
gas and dust to the ISM (about half the total recycled 
gas m aterial) driving further galactic evolution (jMaederl 
Il99l) . 

The mass-loss mechanisms still remain poorly under- 
stood even tho ugh efforts at modeling are decades old 
(|Salpeterl [l974). Thermal pulses enrich the star's outer 
layers by extending a convective zone downward to the 
He s hell and allowing carbon to b e mixed to the sur - 
face (|Schwarzschild fc Harm! 119651 : iHerwig et all l2007f ). 
Material then flows from the stellar surface due to an 
unknown mass-loss mechanism. Once the gas has cooled 
enough to allow the formation of molecules and conden- 
sation of dust grains, the impeding radiation pressure 
from the central source on these grains is sufficient to ac- 
celerate t he grains and drive a dust-driven wind by grain- 
gas d r ag dHovle fc Wickramasinghell962HFleischer et alJ 
119951 iHofner et al.l Il995f ) . The complex interplay be- 
tween the internal thermal pulses, the dust-driven wind, 
and envelope pulsations provides the theoretical frame- 
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work for the observed light curves and mass-loss rates 
(|Whitelock et al.lll99l . 

The efficiency of grain formation and rate of grain 
growth is thus expected to affect the mass-loss history of 
AGB stars. We can constrain dust formation within the 
circumstellar envelopes of AGB stars by analyzing the 
isotopic compositions of certain grains found in primi- 
tive meteorites. Since these meteorite compositions are 
found to differ greatly from solar compositions, they 
must have originated prior to the formation of our solar 
system and are thus dubbed 'presolar.' Detailed stud- 
ies of presolar grain characteristics (such as sizes, crystal 
structures, and composition) can be accomplished after 
extraction from the meteorites (iBernatowicz fc Zinnei 
| 1997tlNittlenl2003l : IClavton k. Nittledl2004t IZinner et al. 
l2006f) . These dust grains are virtually unaltered from 
their formation state in AGB star and thus provide con- 
straints on the physical characteristics of the dust we ob- 
serve in the circumstellar shells of Galactic AGB stars. 
Silcon Carbide (SiC) and carbon graphite dust grains ac- 
count for the majority of discovered presolar grains and 
can be attributed to form ation in Galactic carbon stars 
(|Bernatowicz et aLll2006D . Due to the aggressive and de- 
structive nature of analyzing the meteorites, all but the 
most refractory grains become dissolved. Only recently 
have presolar silicates been discovered meteorites a nd in- 
terplanetary dust particles (jMessenger et al.ll2005[ K 

Initially, the mass loss occurs from the outer layers 
of the star, and hence the chemical composition of the 
recycled material reflects the original ISM composition 
from which the star was formed. As mass loss progresses, 
dredge-up processes enrich the stellar ejecta with nuclear 
fusion p roducts (C, O), s- process elements (Ba, Pb) , and 
lithium (|Drav et aLll2003r ). The nuclear fusion p roduces 
more carbon than oxygen (jLattanzio et alJl2000t l. conse- 
quently causing the C/O ratio of the stellar atmosphere 
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Fig. 1. — Our sample spectra of carbon-rich AGB stars observed in the LMC ranging from 5-36 fira. The spectra have been ordered by 
their [6.4] — [9.3] color from blue to red, normalized by the maximum value of their fitted blackbodies, and offset for clarity. The dashed 
lines indicate the positions of the C2H2 absorption features at 7.5 and 13.8 /im, the SiC emission features at 11.3 /im, and the MgS emission 
feature at ~26 /im. 



to increase from dredge-up of material produced in the 
H- and He -burning layers. For a solar metallicity en- 
vironment, IStraniero et al.l (|1997f ) finds that AGB stars 
above ~1 .5 Mq will eventually become carbon stars un- 
less they exhibit the hot bottom burning phenomenon o f 
more massive stars (M > 5 M Q vr^ lPols et al.ll2001[) . 
However, those with a lower metallicity initially contain 
less oxygen in their atmospheres than their metal-rich 
cousins and, therefore, will more easily become carbon 
stars. Thus, as metallicity decreases, the lower mass limit 



at which stars will evolve into carbon stars also d ecreases 
(|Renzini fc VolilH98ll : IVassiliadis fc Woodlll993h . In ad- 
dition, metallicity may have a profound effect on grain 
growth, the number density of grains, and the grain size 
distribution, thus altering the influence of the radiation 
pressure from the central star. The chemical composi- 
tion of material in the wind will also be affected by the 
metallicity. 

The Magellanic Clouds provide excellent laboratories 
to study the effects of metallicity on the formation of 
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TABLE 1 
Observed LMC Targets 
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Note. — Observed LMC targets sorted alphabetically by the adopted name. RA and Dec coordinates arc given in ,12000. 
Photometry at J, H, and K s are from 2MASS and the A band (8.3 fim) is from the MSX survey. We also list the Astronomical 
Observation Request (AOR) ID along with date of the observation and integration times for the Short-Low (SL) and Long-Low 
(LL) orders in terms of seconds. 
a Skrutskic ct al ct al 



dust. Since the distances to these galaxies and their 
metallicities are well-known, we can derive absolute dust 
mass- loss rates from observed spectral energy distribu- 
tions (|Le Bertrdll997l : lGroenewegen et al.ll2007j ). On av- 
erage, the Large Magellanic Cloud's (LMC) metallicity 
is about half that of the Sun while the metallicity of 
the Small Magella nic Cloud (SMC) is approximately 0.2 
comp ared to solar (Rus sell fc Be sscll 19891 iBertelli et all 
1994). While most Galactic AGB stars appear to show 
an oxygen-rich chemistry, the Magellanic Clouds contain 
a larger fract ion of carbon stars , in line with their lower 
metallicities (jBlanco et al. 198(|). 

Here we present data fro m the Infrared Spectro- 
graph (IRSjHouck ct al. 2004) aboard the Spitzer Space 
Telescope (|Werner et all I2004T ) on 19 carbon-rich AGB 
stars that we co mbine with previously published car- 
bon star samples ijSloan et al.ll206H IZiilstra et al.ll2006t 
lLagadec et al.ll2007t hereafter SZLV The combination of 
these large samples of Galactic, LMC, and SMC spectra 
provides us with a tool to analyze the effects of metallic- 
ity on the composition of the dust and gas and allows us 
to better understand the dust condensation sequence, the 
efficiency of grain formation, and ultimately the mass- 
loss rates. 

2. OBSERVATIONS 

2.1. Sample selection and selection effects 

We present Spitzer IRS data of 19 carbon-rich AGB 
stars in the LMC ( Figure [Tj Table [J) . Based on the color 
classification from lEgan et alj (|2001f ). these objects were 
originally observed for a survey of the dust condensation 
sequence in oxygen-rich AGB stars, but their IRS spec- 
tral signature clearly identify them as being carbon rich. 
This sample covers the same range of H — K s vs K s — A 
color-color space as those LMC objects published by 
IZiilstra et al.l (|2006l ) thereby nearly doubling the avail- 
able data points. The IRS spectrum of IRAS 04496-6958 



was previously published by S peck et al.l (2006). 

In order to analyze the effects of metallicity, we com- 
pare the combined LMC sample to C-rich AGB stars 
observed in the SMC and Milky Way. The SMC stars 
included here were observed u sing the IRS and pre- 
sented in two separat e papers by ISloan et "aLl (|2006l ) and 
lLagadec et al.l (|2007f ). We used Galactic data observed 
with the Infrared Space Observatory (ISO) and pub- 
lish ed in the Short- Wa velength Spectrometer (SWS) At- 
las (jSloan et al.ll2003f ). Thus, we present carbon stars 
from three different metallicity environments: 29 objects 
residing in the SMC, 44 from the LMC, and 34 Galactic 
carbon stars. 

The SWS instrument on ISO observed stars previously 
known to be carbon stars in the Milky Way with a wide 
range of colors. Since the threshold mass to form a car- 
bon star is high in single star evolution in the Galaxy, 
not all of these are 'normal' AGB stars. In particular, 
the blu e stars are more pecu liar high-mass stars or bi- 
naries son fc Richenfl975lh Whereas more than half 
of the Galactic sample contains blue objects expected 
to be naked carbon stars, the LMC and SMC samples 
are much redder on average. This is a known selec- 
tion effect since the Magellanic Cloud surveys were se- 
lected from the mid-IR c olor classification scheme using 
2MASS and MSX colors (lEgan et al.H200l . Since MSX 
is only sensitive towards objects with a strong 8 /zm ex- 
cess, it was expected that most of the observed objects 
would have high mass-loss rates and relatively red spec- 
trum. Indeed, for most of the Magellanic Cloud objects, 
we find gas mass-loss rates of 10 -6 M Q yr _1 or higher 
while Galactic AGB stars have a wider range of 10 -8 ' 5 
to 10~ 4 - 5 M Q yr" 1 (Section EES}. The Galactic sam- 
ple contains the widest range of colors spanning between 
0.3 < H - K s < 4.5 while the LMC and SMC objects 
range between 1 < H - K s < 3.2 and 0.9 < H - K s < 2.1 
respectively. On average, the carbon stars observed in 
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the LMC are redder than those observed in the SMC. 

2.2. Spectral features 
2.2.1. C2H2 absorption 

For the majority of stars in the three samples, the most 
prominent molecular absorption features occur at ^7.5 
and ~14 /im and are normally associated with acetylene 
(C2H2) absorption. However, HCN and CS absorption 
are also know n to provide contributions to the dep th of 
these features (|Goebel et a l.lll98l UAoki et al.lll998 ). Us- 
ing the 3.5 /im absorption feature. Ivan Loon et alj (|2006l ) 
show that HCN and CS are compar atively weak for LMC 
stars while iMatsuura et al.l ()2006l ) found no clear evi- 
dence of HCN among their LMC sample. In particular, 
they interpret the low HCN abundance in the Magellanic 
Clouds as being due to initially low nitrogen abundances 
and effi cient carbo n prod uctio n in metal-poor e n viron - 
ments. lAoki et al.1 (|1999D and ICernicharo et alj (|1999h 
only detect HCN in 1/3 of their Galactic carbon sample. 
Thus, we can confidently identify these two absorptions 
as prominently due to C2H2. 

The double peaked absorption band centered at 7.5 /im 
arises from the C2H2 v\ + v\ transitions in the P and R 
branches. The Q branch band of the C2H2 fundamental 
V5 bending mode creates a sharp absorption feature at 
13.7 /im. Colder gas in the extended envelope (possi- 
bly above the dust-forming region) produces the sharp 
Q-branch feature while the broad P and R branches 
originate from warmer gas deeper into the atm osphere 
(jJorgensen et al.ll2000f ). IMatsuura et all (|2006h explain 
the broad wings extending on either side of the 13.7 /im 
feature as being due to the P and R branches. In the most 
extreme cases, these wings can extend between ~11 and 
~17 /um (jSpeck et alJl2006ft . Throughout this paper, we 
refer to the Q branch as the "13.7 (im feature" and the 
P+R branches in the same region as the "14 /im feature" . 

2.2.2. 10-micron absorption 

A significant number of stars exhibit a relatively weak 
absorption band spanning ^8.5 to ~11 /im with a peak 
depth around 9.4 /im. We mainly observe this feature 
in the bluest objects. While the carrier of this feature is 
still in doubt, our observations support a molecule that 
originates near the photosphere of the star and would be 
masked by the infi ll ing of dust emission in redder o bjects. 

iJoTgensen etaD (|2000t ) and lZijlstra et all (|2006f ) favor 
an identification of C3; however, we do not observe a 
corresponding absorption between 6.6 and 7 /im that is 
also predicted. Furthermore, we cannot use the deep ab- 
sorption feature at 5 /im to confirm the presence of this 
molecule since this band also contains significant contri- 
butions from CO. Laboratory measurements are neces- 
sary in order to firmly identify C3 as the carrier of the 
10-micron absorption. 

2.2.3. SiC emission 

The majority of the carbon stars in the Magellanic 
Clouds and the Galaxy show the presence of an emission 
featu re around 11.3 /im gener ally attributed to SiC emis- 
sion (jTreffers fc Cohenlll974| ). While the central wave- 
length of the SiC resonance as measured in the laboratory 
intrinsically lies at 11.3 /im (|Gilralll97lj ). the central peak 
in astronomical data varies between 11.15 and 11.7 /im 



for Galactic sources. ISpeck et al.l (|2005f ) used grains of 
different sizes to interpret this shift: larger grains pro- 
vide a cool absorbing component around 10.8 /im while 
a warmer emitting component towards the red arises 
from smaller grains. Thus, different grain-size distri- 
butions will theoretically provide different contributions 
from emitting and absorbing components, causing the 
overall central wavelength of the SiC feature to shift. 

2.2.4. MgS emission 

In slightly less than half of all the objects presented 
here, we observe a broad emission feature beginning blue- 
ward of 22 /im and extending beyond the IRS spec- 
tral range. This feature is commonly seen in carbon- 
rich AGB and post-AGB stars as well as planetary neb- 
ulae and consists o f two subpeaks at 26 and 33 /im 
([Forrest et al.1 Il981h both of which are attribute d to 
magnesium sulfide (MgS; iGoebel & Molielevl ll985h By 
calcula ting absorption cr oss sections for various grain 
shapes. iHonv et all (|2002D modeled the MgS spectral fea- 
ture of Galactic carbon stars and showed that the sub- 
peaks are due to different distributions of grains shapes 
and that the peak wavelength of the feature shifts with 
temperature. 

3. ANALYSIS 

3.1. Manchester method 

Since the dominant component of the dust in 
these sources consists o f featureless amorphous carbon 
(jMartin fc Rogersl [l987). it is useful to develop an ob- 
servable quantity that represents the dust opacity. Pre- 
vious analyses of car bon star data ha ve utilized the 
"Manchester method" (jSloan et al.ll2006h in order to de- 
fine the amorphous carbon dust continuum amid the 
plentiful spectral emission and absorption features. This 
method defines four narrow bands centered at 6.4, 9.3, 
16.5, and 21.5 /im that are used to determine a set of 
two color temperatures of [6.4] — [9.3] and [16.5] — [21.5]. 
These colors can then be used as proxies for the optical 
depth of the surrounding dust shell (r oc [6.4] — [9.3]) and 
for estimating the dust temperature, respectively. 




Wavelength (/im) 



Fig. 2. — The spectra of U Cam showing the effect of the 10 fim 
absorption feature on the [9.3] Manchester band represented as 
the cross-hatched vertical bar. The 10 fim band gives a flux value 
significantly below the probable continuum level, which is approxi- 
mated by a modified black body. The resulting straight line contin- 
uum fit (solid line) differs significantly from the modified blackbody 
fit (dashed line) underneath the SiC emission. 
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Fig. 3. — Example spectra with error bars of three sources, one 
from each environment - SMC, LMC, and the Milky Way from 
top to bottom, respectively. The dotted line indicates the chi- 
square fits of the underlying carbon-rich dust continuum using a 
single-temperature modified blackbody. The straight line segments 
represent the estimated continuums using the method employed by 
SZL. We find the feature strengths by dividing the derived black- 
body fits. 

A few complications arise when trying to perform this 
method on a sample of objects containing a wide range of 
dust temperatures. In particular, the Manchester 9.3 /im 
band fall directly within the 10 /mi absorption feature, 
thus underestimating the dust continuum flux at this 
wavelength for stars with this feature (Figure [2]). How- 
ever, this only affects a small number of the bluest ob- 
jects providing a mean shift of 0.07 for the [6.4] — [9.3] 
color. Thus, we still rely on the [6.4] — [9.3] color to 
provide us with a qualitative analysis of the dust opaci- 
ties. Furthermore, the 21.5 /im band, used to define the 
blackbody continuum, can acquire increased flux from 
the adjoining broad MgS feature for significant emission 
strengths, causing the continuum emission to be overes- 
timated at these wavelengths. Prominent MgS emission 
can occur shortwa rd of 20 /im for low dust temperatures 
dHonv et al.ll2002Th 

Additionally the Manchester method adopts straight- 
line continuum fits in order to measure the strengths of 
the C2H2 absorption bands and the SiC emission feature. 
They defined wavelength segments of supp osed contin- 
uum that flank the features of interest (see iSloan et"aT1 
2006). However, since the blue side of the continuum 
used for the SiC feature ranges from 9.5 to 10.1 /im, 
this method will tend to underestimate the underlying 
continuum and overestimate the flux in the SiC band 
(Figure [2]). Similarly, the varying strength of the C2H2 



feature affects determination of the red side of the SiC 
feature. In the most extreme cases, C2H2 absorption can 
reach as low as ~11 /im causing considerable suppres- 
sion of the observed SiC emission. The same variation 
in broadness of the C2H2 bands also makes it difficult to 
set a consistent wavelength range at 7.5 and 14 /tm that 
allows for proper isolation of the P, Q, and R branches 
of C2H2. For this reason, previous studies only focused 
on the narrow Q branch of the 13.7 /im feature and dis- 
regarded the flanking P and R branches. While this pa- 
per strives to improve on quantitative determinations of 
the various spectral features, we stress that our improve- 
ments described in the following section do not change 
the majority of the qualitative conclusions found by SZL 
using the Manchester method. 

3.2. Modified blackbody fits 

SZL employed the Manchester [16.5] — [21.5] color to de- 
rive a blackbody color temperature and thereby extrap- 
olate the continuum underneath the MgS feature. How- 
ever, the aforementioned problems associated with this 
type of continuum estimation can lead to inaccurate de- 
terminations of dust properties such as the temperature, 
mass-loss rate, and total dust mass. We improve on this 
method by fitting a modified blackbody to each carbon 
star's spect r um. E xpanding on the method described by 
iHonvet all (|2002'1. we utilize a x fitting routine to find 
the best single temperature modified blackbody to model 
the entire spectral continuum, 

F{\) =A x B{X,T) x \~p, (1) 

where F(X) is the flux density of the continuum at some 
wavelength A, B(X, T) is the Planck function with a tem- 
perature T, p is the dust emissivity index, and A is a scal- 
ing factor. We minimize the x 2 values by varying the free 
parameters (T, p, and A) from Equation [1] over selected 
continuum points. The p-parameter represents how effi- 
ciently the dust grains radiate at wavelengths larger than 
the grain size. To some extent, varying this parameter 
allows us to incorporate effects from optical depth and 
temperature gradients. Crystalline material has p ~ 2, 
amorphous material has an emissivity index between 1 
and 2, and layered materials have p ~ 1. Temperature 




oL , , , 1 , , , , 1 , , , , 1 , , , J 

0.0 0.5 1.0 1.5 2.0 

[6.4]-[9.3] 

Fig. 4. — Contours of constant p-parameter values for modified 
blackbodies with temperatures ranging between 300 and 2500 K 
with —1 < p < 1.4 with a 0.2 increment. The symbols represent 
data from fitted blackbodies for the samples of Galactic, LMC, and 
SMC objects. Error bars for [6.4] — [9.3] are on the order of the 
symbol size. 
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Fig. 5.— Correlation between 2MASS J - K and [6.4] - [9.3] 
colors for Galactic objects. We have overplotted the best-fit linear 
relation. 

gradients within the dust shell lower the value of p by 
broadening the spectral energy distribution. Similarly, 
an optically thick dust shell further reduce the value of 
P- 

For a source embedded in a circumstellar shell, one ex- 
pects the continuum emission to follow the general form 
of 



F V (X) = J v (A,T*)e- T - +B v (\,T os )(l 



(2) 



where T+) is the stellar emission, d is the distance to 
the source, B u (X,Tcs) is the blackbody emission of the 
circumstellar dust shell, r„ is the optical depth of the dust 
shell, and and Tcs are the prevailing temperatures of 
the stellar photosphere and envelope, respectively. In 
cases where the optical depth is large at infrared wave- 
lengths, the contribution from the stellar photosphere 
in Equation [3 is negligible, and the modified blackbody 
fit (Equation [1]) gives an approximation for the average 
temperature of the dust shell. However when the op- 
tical depth is small, the photospheric emission becomes 
important, and both and Tcs contribute to the tem- 
perature in Equation[T] For this case, the temperature T 
from the modified blackbody fit should be treated only 
as a fit parameter. 

The blackbody continuum fits allow us to improve on 
straight line continua by determining more precisely the 
wavelength range over which spectral features are present 
(Figure [3]). In order to accommodate for the varying 
widths of the features, the continuum points were se- 
lected for each spectrum individually. We chose contin- 
uum points redward of the 7.5 /im C2H2 band at; be- 
tween the 7.5 /Ira C2H2 band and SiC emission feature; 
and between the 14 /im C2H2 and MgS feature, taking 
care that none of the feature wings were included in the 
continuum selection. When fitting ISO spectra, the in- 
clusion of wavelengths shortward of 5.2 firm that are un- 
available in IRS spectra alters the best-fit temperature 
up to 10% and slightly shifts p towards negative val- 
ues. In the majority of these cases, though, only minor 
differences arise in the temperature and p-value when 
considering the 2.3 to 5.2 /im region to constrain the 
continuum. Moreover, changes to the fitted continuum 
above and below the features of interest were impercep- 
tible; therefore, we do not introduce a bias in the feature 
strengths of Galactic objects. 

In addition, we considered 2MASS J, H, and K s points 



(|Skrutskie et al.l l2006f ) to constrain the shorter wave- 
lengths of the blackbody points. However, due to the 
intrinsic variability of these objects, the 2MASS data 
points were only used to guide the slope of continuum fit 
rather than as rigid anchors. Additionally, different dust 
temperatures dominate the emission between the near- 
IR 2MASS photometry and the mid-IR spectra, but we 
partially account for this effect with the inclusion of the 
dust emissivity index, p. Also, due to the relative dis- 
tance and faintness of objects in the Magellanic Clouds 
as compared to objects in the Galaxy, spectral and pho- 
tometric data for LMC and SMC objects contain larger 
uncertainties than Galactic objects. In the entire sam- 
ple, 28 objects have a J-band magnitude greater than 
the lO-er point-source detection limit of 15.8 mag, with 
23 of these objects residing in the LMC. This number 
improves by half for the 10-tr H-band limit of 15.1 mag. 
Thus, using the 2MASS points as absolute constraints 
for the single-temperature blackbody fits was not recom- 
mended in every case. 

Figure 2] shows the result of Equation [1] for modified 
blackbodies ranging from temperatures of 300 to 2500 K 
and p-values of —1.0 to 1.4. Continuous lines show how 
the temperature changes for varying [6.4] — [9.3] with 
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Fig. 7. — From left to right, the top panels plot the equivalent widths for the C2H2 P and R branches at 7.5 fim as a function of 
[6.4] — [9.3] color (dust mass-loss rate) and gas mass-loss rates, respectively. The bottom panels show the equivalent widths for C2H2 P 
and R branches at 14 /an as a function of [6.4] — [9.3] and gas mass-loss rates. 



constant p-values. We overplot the derived p and T 
values against the measured [6.4] — [9.3] values for the 
SMC, LMC, and Galactic samples. The mean p- values 
for the entire sample is 0.13±0.26, and for the individual 
samples are 0.27 ± 0.2 (SMC), 0.18 ± 0.22 (LMC), and 
—0.05 ± 0.25 (Galactic). The variation in average values 
between the three environments may indicate a metallic- 
ity effect; however, since the p-parameter represents the 
optical properties of the grain which evolve over the life- 
time of these objects, this may reflect the biases of each 
sample. For instance, the objects within the Milky Way 
contain the widest range of colors with the majority re- 
siding towards the blue, while the SMC sample contains 
the smallest color range and consists of much redder ob- 
jects. 

Due to the poor quality of some spectra, mis-pointings, 
and/or background contamination, we were unable to 
obtain blackbody fits for a total 11 objects from the 
various samples: MSX LMC 78 2 from our sample ; 
NGC 1978 MIR1, NGC 1978 IR4 (|Ziilstra et al.ll20M k 
NGC 419 L E 18, NGC 419 LE 2 7, NGC 419 LE 35, and 
ISO 00549 (|Lagadec et a l.l l200l: and IRAS 18 430-0237. 
S Lyr, SS73 38, and SZ Sgr (|Sloan etlj1l2"003h . 

Table 2 lists the feature strengths determined using our 
modified blackbody fits. We measure equivalent widths 
of the 7.3 /im band, 13.7 /im C2H2 Q branch, and the 
P+R branches at 14 /im. Strengths of the dust emis- 
sion features were determined by taking the ratio of the 
integrated flux to the integrated continuum. The cen- 
tral wavelength for the SiC feature was also found. We 



do not include features where low signal-to-noise ratios 
cause uncertainties comparable to the strength. 

3.3. Mass-loss Rates 

Based on du st radiative transfer models, 

iGroenewegen et ail (|2007h determined dust and gas 
mass-loss ra tes for 60 carbo n stars in the Magellanic 
Clouds while [Le Bertrd (|1997l ) modeled the circumstellar 
dust shells of 23 Galactic carbon stars. These works 
derived Mdust directly from the infrared SEDs and NIR 
photometry. Gas mass-loss rates were subsequently 
derived by scaling by an assumed dust-to-gas ratio (tp) 
based on the metallicity. Even though this is a common 
assumption when deriving M gas , there are no studies 
that conclusively show that the total mass-loss rat e 
explicitly depends on metallicity (|Ziilstra et al.1 [1996T) . 
The dust-to-gas ratio has been shown to depend approx- 
imately linearly on initial metallicity while mass-los s 
rate shows only a weak dependence (jvan Loonl[2000Tl . 
However, in order to proceed with our analysis of the 
dependencies of dust and molecular abundances with 
pressures and densities, we must make some assumption 
about the gas mass-loss rates. 

When available, we use m ass-loss rates pub lished by 
IGroenewegen et all (|2007l) or lLe Bertri (|1997l ). In their 
derivation of M gas , IGroenewegen et all ([2007D used the 
same dust-to-gas ratio of ip = 0.005 for all objects 
within the LMC and SMC. Instead, we have adopted 
ip = 0.002 for SMC objects, corresponding to a metal- 
licity of Z/Zq = i. This assumes that ip scales with 
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Fig. 8. — (a): Plot of the equivalent with for the 13.7 (im C2H2 Q branch as a function of the [6.4] — [9.3] color (dust mass-loss rate), 
(b): The C2H2 equivalent width for the 13.7 um Q branch as a function of the derived gas mass-loss rate. 



metallicity and ipMW = 0.01. 

For the obje ct s w ithin Magellanic Clouds, 
iGroenewegen et "all ([2007) show a correlation be- 
tween M dust and the [6.4] - [9.3] color (Figure Ek) . The 
gas mass-l oss rate for the LMC objects can then be 
derived by (|van Loonll2007ft : 

log(M gaSiLMC ) = 2.55([6.4] - [9.3]) ' 5 - 7.6. (3) 

For the SMC sample, we modify Equation[3]using a value 
of Vsmc = 0.002, 

\og{M gas ,sMc) = 2.55([6.4] - [9.3]) ' 5 - 7.2. (4) 

All mass- loss rates are in M Q yr -1 . 

Since these relations were derived only for objects be- 
tween 0.2 < [6.4] — [9.3] < 1.5, the mass-loss rates of 
objects with [6.4] — [9.3] < 0.2 are not well constrained, 
affecting 11 Galactic objects. Instead, we determined gas 
mass-loss f or all Galactic objects using the relation from 
He Bertrd (fl99l . 

\og(M gasM w) = -6.0/(J - K - 0.2) - 4.0, (5) 

assuming a dust-to-gas ratio of 0.01. We find a linear 
correlation between the J — K and [6.4] — [9.3] colors 
for Galactic objects ( Figure Et. There fore, we can over- 
plot the trend from iLe Bertra (|1997f ) in terms of the 
[6.4] — [9.3] color. The top panel of Figure [6] shows that 
the dust mass-loss rates of the Galactic objects derived 
from 2MASS J — K indeed overlap with those derived 
from the [6.4] — [9.3] color. Once we apply metallicity 
corrections for determining the gas mass-loss rates, ob- 
jects in the separate samples no longer fall along a single 
trend (Figure [6p). Table [3] lists the dust and gas mass- 
loss rates for each object. 

4. RESULTS 
4.1. Acetylene molecular absorption 
4.1.1. P & R branches 

We plot the equivalent widths of the C2H2 bands at 7.5 
and 14 /im as a function of the [6.4] — [9.3] color and M gas 
for the Galactic, LMC, and SMC samples in Figure \7\ 
By plotting against the [6.4] — [9.3] color, we are com- 
paring sources with the same optical depth. Assuming 
a constant expansion velocity independent of metallic- 
ity, the [6.4] — [9.3] color also directly correlates with the 



dust mass-loss r a te ass uming (jGroenewegen et al"1l2007t ). 
iMatsuura et alj (|2002f) show that increasing the contri- 
bution from dust emission can affect the molecular band 
strength by filling in the C 0H2 absorption and dec reas- 
ing the equival e nt wi dth. Ma tsuura et al.l (|2005f) and 
Ivan Loon et al.l (|2006l ) show evidence of this "dust di- 
lution" for C2H2 bands occurring at 3.1 and 3.8 /im in 
the Magellanic Clouds. 

Analyzing the Galactic sample, we find that the 7.5 /im 
band shows a behavior simil ar to the 3.8 /im C2H2 bands 
from Matsuura ct al. (2005). Figures [7K and [7b show an 
initial increase in equivalent width with optical depth 
and mass-loss rate, indicating a rise in level popula- 
tions from increasing pressure and density. Circumstel- 
lar emission subsequently fills in the absorption band at 
Mgas ~ 10 -6 5 causing a decline in the equivalent width. 
The 7.5 /im absorption feature reflects this trend better 
than the 14 /im band since the SiC emission feature tends 
to overlap the 14 /im P branch. 

The P and R branch acetylene bands also show a de- 
pendence on metallicity. For a given [6.4] — [9.3] color, 
objects from the SMC sample show slightly larger equiva- 
lent widths than those from the LMC. This is even more 
pronounced when compared to Galactic sources at 7.5 
and 14 /im albeit with considerable scatter (Figures [7^ 
and [7J:). However, plotted against gas mass-loss rate, 
those objects from the LMC and SMC occupy the same 
space (Figures [7b and[7ji). Our assumptions when deriv- 
ing gas mass-loss rates effectively discard the metallicity 
de pendence between th e LMC and SMC objects found 
bv lLagadec et al.l (|2007l ) when using different and equally 
valid assumptions. Thus, we must be aware that the as- 
sumptions we make affect qualitative results when deal- 
ing with the gas mass-loss rates. However, the LMC and 
SMC objects show significantly higher equivalent widths 
than their Galactic counterparts with the same M gas as 
found by SZL. These works suggest that efficient dredge- 
ups for low metallicity C-stars help to increase the C/O 
ratio and, therefore, the C2H2 equivalent width for these 
bands. Additionally, for a given M gas , objects in the 
Magellanic Clouds will contain less dust than Galactic 
objects causing effects from dust dilution to be less pro- 
nounced for the SMC and LMC samples. 

4.1.2. 13.7 urn Q branch 
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Fig. 9. — (a): Correlation between the 14 fim equivalent width and SiC central wavelength. Error bars for Galactic objects are smaller 
than symbols. Objects with increasingly large C2H2 band absorptions tend to have XsiC shifted further towards the blue while objects 
with redder \siC have smaller C2H2 equivalent w idths, (b): The SiC apparent central wavelength as a function of the [6.4] — [9.3] color 
similar to Figures 9 & 10 from lSloan et all ||2006T) . The dashed lines are linear fits to trends described in the text. The horizontal dotted 
line indicates the SiC central wavelength observed in laboratory data. 



Figure [5] shows the equivalent width of the Q branch 
band at 13.7 /mi as a function of the [6.4] — [9.3] color 
and Mgas- For Galactic objects, we observe a trend 
of the 13.7 /mi Q branch equivalent width with mass- 
loss rate. The equivalent width initially decreases until 
[6.4] - [9.3] = 0.5 and M gas = 1(T 6 M yr" 1 then very 
slightly increases. Van Loon et al. (2006) observe similar 
trends in absorption bands at 3.1, 3.57, 3.7, and 3.8 /mi 
and explain them as being due to a competition between 
the effect of dust dilution and the effect of higher densi- 
ties from an increasing mass-loss rate. The initial decline 
towards red colors in the Galactic sample suggests that 
circumstellar dust emission fills in the absorption fea- 
ture. As gas mass-loss rate and pressure increase, densi- 
ties become sufficiently high to overcome effects from the 
veiling at around 10 -6 M Q yr -1 and effectively increase 
the 13.7 /mi equivalent width. This would suggest that 
the 13.7 /mi C2H2 band has a stronger dependence on 
the gas density than the dust column density at higher 

Mgas- 

We observe a similar trend of initially decreasing 
13.7 /zm equivalent width for SMC objects in Figure 
at [6.4] — [9.3] < 0.5, which overlaps the Galactic trend. 
However, when plotted against M gas , these particular 
objects no longer correspond to the observed Galactic 
trend, supporting the interpretation that dust dilution 
has a prominent influence on the 13.7 /mi absorption 
at lower mass-loss rates. Such a trend is not observed 
for the LMC sample since only one object (GRRV 38) 
has a value of [6.4] — [9.3] less than 0.5. At M gas > 
10~ 6 M Q yr -1 , though, the slightly increasing trend of 
the Galactic objects is generally matched by both the 
LMC and SMC objects. This further suggests that effects 
from increasing densities dominate over effects from dust 
dilution at higher gas mass-loss rates for this absorption 
band. 

Figure 5 of iMatsuura et all (|2006D show the same in- 
crease in the 13.7 /mi equivalent width with [6.4] — [9.3] 
color and interpret this as evidence for a lack of dust 
dilution. However, except for one star, their entire sam- 
ple consists of objects redder than [6.4] — [9.3] = 0.5 



with M gas ^ 10~ 55 Mq yr -1 ; therefore, they were 
unable to observe the initial decrease in the equiva- 
lent width. Although we observe the more obvious ef- 
fects of dust dilution at lower optical depths and lower 
M, w e still support the conclusions of IMatsuura et all 
(2006) that the molecules exist throughout the envelope. 
Since the equivalent widths do not continue to decline 
at higher mass-loss rates where photospheric molecules 
would be veiled by the overlying dust, then the absorbing 
molecules most likely exist throughout the circumstellar 
envelope. 

4.2. Dust emission 
4.2.1. SiC central wavelength 

The central wavelength of the SiC feature (Xsic) nas 
been us ed to infer dust pr operties and circumstellar ge- 
ometry ([Speck et all 12005). but is very sensitive to the 
placement of the continuum. Despite the improved con- 
tinuum determination, spectra with prominent acetylene 
bands at 14 /mi continue to affect the central wavelength 
by removing emission from the SiC feature (Figure [5^) . 
Objects with larger XsiC tend to have little C2H2 absorp- 
tion while objects with significant C2H2 absorption have 
K<n c occurring towards the blue. IZiilstra et alj (|2006f ) 
and lLagadec et all (|2007t ) raise the possibility that since 
the C2H2 absorption overlaps the SiC emission, the red 
wing of the SiC feature is suppressed and causes the cen- 
tral wavelength to shift towards shorter wavelengths (see 
also Figure fTUk) . 

C2H2 suppression can explain the initially decreasing 
trend for Galactic stars when p l otting [6.4] — [9.3] against 
XsiC reported by ISloan et "all (|2006r) and displayed in 
Figure [Sb. Most of the blue Galactic objects have little 
C2H2 absorption; therefore, the SiC feature has not been 
altered on the long-wavelength side. At the same time, 
absorption at ~10 /mi lowers emission from the blue side 
of the SiC feature, as shown in the Figure [TUb . This 
has the opposite effect of the C2H2 absorption, causing 
Xsic to be artificially shifted to longer wavelengths as 
the temperature increases ([6.4] — [9.3] color decreases). 
Thus, for the bluest objects, we observe Xsic to be redder 



10 



Leisenring et al. 



S -1 




11 12 13 14 

Wavelength (/itm) 



10 11 12 

Wavelength (fim) 



Fig. 10. — (a): The dark line shows the normalized continuum subtracted spectrum of LI-LMC 1817. Ovcrplotted in gray is the spectrum 
of MSX LMC 1298. C2H2 absorption including P, Q, and R branches are labeled as well as the SiC emission feature. The vertical dotted 
lines indicate the centroid of the SiC feature as it is affected by C2H2 absorption, (b): The spectrum of DO 16793 is plotted in black while 
RY Dra is in gray. The 10 /^m absorption and SiC emission features are labeled, and the vertical dashed lines show the shift of \sic when 
affected by the absorption band. 



than the lab data. 

For [6.4] — [9.3] > 0.4, Galactic carbon stars show in- 
creasing Xsic with increasing values of [6.4] — [9.3]. SMC 
stars appear to follow the same general trend as the 
Galactic sample over the observed color range, albeit a 
little less pronounced (Figure [9Jd). However, the sample 
of SMC objects contains fewer stars overall and spans a 
much smaller color range. LMC stars, on the other hand, 
show a fairly consistent central wavelength between 11.1 
and 11.2 fj,m for [6.4] - [9.3] = 0.5 to 1.5. These ob- 
jects contain some of the most substantial 14 /zm C2H2 
bands which dramatically affects the observed SiC cen- 
tral wavelength. 

The different trends of the SiC central wavelength 
observed for the LMC, SMC, and Galactic samples at 
[6.4] — [9.3] > 0.4 can all be explained by a competi- 
tion between self- absorption of the SiC feature and ef- 
fects from C2H2 absorption. SiC self-absorption shifts 
the feature from 11.3 jum towards the red while the broad 
14 fj,m C2H2 absorption band causes the SiC feature to 
shift towards the blue. Hence, objects containing little 
C2H2 absorption show a steeper slope with increasing op- 
tical depth (ie. Galactic sample). Since the Magellanic 
Cloud samples contain significantly wider C2H2 bands, 
the slope of this trend decreases. 

Two objects, AFGL 341 and IRAS 21489+5201, stand 
out with extremely red central wavelengths at about 
12.1 /jm. These particular objects are among the reddest 
sources, exhibit no apparent 10 jum absorption feature, 
and contain a relatively weak C2H2 absorption at 14 /zm. 
Thus, self-absorption of the SiC feature most likely pro- 
duces the considerable shifts in the central wavelength 
for these two objects. 

While self-absorption of the SiC feature may have a 
considerable effect on trends of the Galactic sources, it 
is difficult to make a conclusive statements due to the 
molecular absorptions flanking the dust feature. We find 
only two objects in our entire sample that exhibit a shift 
in Xsic mainly due to SiC self-absorption. Therefore, 
one must be careful when making claims about the state 
of the SiC optical depth and self-absorption based purely 
on the position of its central wavelength, especially in the 



presence of strong molecular bands. 

4.2.2. SiC feature strength 
For a given [6.4] — [9.3] color (which correlates with 
dust optical depth and Mdust), the line-to-continuum ra- 
tio of the SiC emission feature appears largest for Galac- 
tic objects, but shows no appreciable difference between 
the LMC and SMC objects over their overlapping color 
range of 0.5 < [6.4] - [9.3] < 1.0 (Figure [TTJi). Plotting 
the SiC strength against the gas mass-loss rate provides 
a good indication of the observational bias inherent in 
the Magellanic Cloud samples (Figure [Tib). GRRV 38 
and RAW 960 are the only two sources within the Magel- 
lanic Clouds with M gas < 10~ 6 M Q yr" 1 , while about 20 
Galactic objects span the range between 10 -9 and 10 -6 
Meyi" 1 . 

While trends within the Magellanic Clouds are difficult 
to discern due to the large scatter of the SMC sources 
and small number of LMC sources with [6.4] — [9.3] < 0.5, 
the Galactic sources increase rapidly and pe ak at [6.4] — 
9.3] ~ 0.5, matching the observations by ISloan et ail 
20061) . This same trend also transfers to Figure [Tib 
with a general increase in SiC emission relative to car- 
bon emission from low (10 -9 M Q yr _1 ) to intermediate 
(10~ 6 - 5 Mq yr" 1 ) mass-loss. 

Moving towards redder colors and higher mass-loss 
rates, the SiC strength subsequently decreases. The 
starting point of this decline appears to depend on metal- 
licity. This trend is most pronounced in the Galactic ob- 
jects which initially increases, then plateaus, and finally 
decreases between 10~ 9 to 10 -4 - 8 M yr -1 . The same 
decreasing trends for the LMC and SMC samples occur 
at M gaa ~ 10" 5 - 3 and M gas ~ 10~ 5 M© yr" 1 respec- 
tively and line up well with the Galactic objects. 

The trend of decreasing SiC with increasing color could 
result from self-absorption within the feature due to an 
outer layer of cooler dust rather than a decrease in actual 
SiC grain emission ([Speck et al J 12005ft and may explai n 
the observations of IRAS 04496-6958 (Speck et al.ll2006h . 
However, this trend also correlates with the onset of MgS 
emission and may be more closely linked to the dust con- 
densation sequence (see Sections 14.2.31 and 5). 
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Fig. 11. — Plots of the SiC and MgS dust feature strengths with respect to the [6.4] — [9.3] color (dust mass-loss rate) and the gas 
mass-loss rate. Objects with a measurable MgS emission feature have been marked in the top two panels with a black X. In general, redder 
sources have a weaker SiC emission and correlate with the presence of MgS. Dotted lines in all panels indicate the transition point in each 
metallicity system from objects without any MgS to objects with MgS emission. 



4.2.3. MgS feature strength 

Every carbon star throughout the entire range of color- 
temperatures presented here contains an apparent SiC 
dust feature; however, less than half of the entire sample 
show a MgS band. Only six out of 29 carbon stars in the 
SMC sample exhibit a MgS feature, 27 out of 44 in the 
LMC, and 14 out of 34 for Galactic sources. 

Throughout the entire sample, there appears to be an 
overarching trend for redder objects to contain a stronger 
MgS emission relative to the continuum. All objects with 
a [6.4] — [9.3] color greater than 0.9 (corresponding to a 
Tcont ^$ 500 K) contain a significant MgS dust feature 
while all objects with [6.4] — [9.3] < 0.6 show no emission 
from MgS, with the exception of V Cyg, V CrB, and 
T Dra in the Milky Way. This is consistent with the 
condensation temperature of MgS (~600-300 K), which 
is always lower than the SiC condensation temperature. 

We also find that the MgS condensation temperature 
depends on pressure. Plotting the MgS strength in terms 
of the gas mass-loss rate reveals that the majority of ob- 
jects with an MgS feature lie within a relatively small 
mass-loss range, albeit a bit offset in terms of metallic- 
ity. About half of the Galactic spectra containing an 
MgS feature have continuum dust te mperatures greater 
than t he condensation temperature. iHonv fc Bouwmanl 
(2004) have explained this phenomenon in terms of an 
optically thin detached dust shell due to time-variable 
mass loss. Such objects are identified by a blue spec- 
trum, low mass-loss rate, and weak MgS feature. 



Each metallicity environment appears to contain a 
sharp transition region where MgS emission is observed 
which occurs at different mass-loss rates and color tem- 
peratures (vertical dotted lines in Figure [TTj) . These off- 
sets can be attributed to metallicity differences and can 
be used to constrain the conditions under which we would 
expect to observe MgS grain emission, assuming the MgS 
and amorphous carbon exist in the same volume (ie. no 
detached dust shell). 

These transition points also correspond with the de- 
crease in SiC emission discussed in Section 14.2.21 Fig- 
ures [T3Ji andfTTb show the correlation between SiC emis- 
sion strength and the MgS feature with respect to the 
dust and gas mass-loss rates. Objects containing weak 
MgS features mainly occupy positions near the peaks 
of the SiC trend discussed above. As MgS intensity 
strengthens, SiC declines suggesting that the MgS con- 
denses in the form of a coating on top of the SiC grains, 
thus suppressing the SiC emission (see Section 5). 

5. DUST CONDENSATION SEQUENCE 

In order to decipher the dust condensation sequence 
in carbon stars, we must first constrain various param- 
eters in the envelope such as temperature, gas pressure 
(density), C/O ratio, and abundance (metallic ity). For 
solar metallicity stars, lLodders fc Feglevi (1995) give the 
condensation temperatures of C, SiC, and MgS for vary- 
ing C/O ratio and pres sure (Fig ure fT2"k). At a given 
C/O ratio and pressure. iLoddera (|2006f ) shows that the 
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Fig. 12 , — Condensation temper atures of carbon graphite, SiC, and MgS as a fu nction of t he tot al pressure and C/O ratio (1.05 to 2.0) 
based on Loddcrs k Fcglcy (1995). Data for the LMC and SMC are scaled from Loddcrs (2006). The dash-dotted line shows a model 
atmo sphere for a Carbon AGB star with the following parameters: M* = 1 Mq, L* = 10 4 Lq, T e g = 2400 K, and C/O = 1.7 l|Woitkel 
120061 1. The filled circles correspond to a radius of 1.1 il*, 1.5 i?*, 3 i?*, and 6 R* from top to bottom, respectively. The curves in panels b 
and c have been shifted to compensate for metallicity effects. 



condensation temperatures will decrease with decreasing 
metallicity; however, this effect is more pronounced for 
carbon than for SiC or MgS. Figures [12b and [12b plot 
the estimated shift in the condensation temperatures for 
metallicities corresponding to the LMC and SMC, re- 
spectively. We also compare th e radial structure of a stel- 
lar model from IWoitkd (|2006f ) to show the approximate 
regions a nd sequence of various dust condensates. Fur- 
thermore, Jorgcns en et alj (|1992f ) show that for a given 
mass, decreasing the effective temperature or the C/O 
ratio will generally shift this curve downward and to the 
left. A reduction in the metallicity similarly reduces the 
pressure at a given radius, but the temperature structure 
will remain fairly constant. 

Observationally, prominent MgS does not appear until 
the mass-loss rates for Galactic, LMC, and SMC objects 
become greater than lCT 58 , 1CT 54 , and 10~ 51 M© yr" 1 , 
respectively (Fi gure [TTB) . We , therefore, support the 
conclusions frorn lLoddersI (|2006t ) that for a given temper- 
ature, decreasing the metallicity also leads to an increase 
in the critical density at which MgS condenses. 

Figure [13] is a schematic of the trends we expect to see 
for three different dust condensation sequences to which 
this method is sensitive. These can be differentiated by 
analyzing the change in the SiC-to-C ratio with M gaa and 
taking into account the formation of MgS (Figure [TTJd): 

I. C — > SiC — > MgS: In general, carbon stars will start 
with a low C/O ratio and a low mass-loss rate. As 
the wind expands away from the star and the tem- 
perature decreases, carbon grains condense first. 
As the temperature further decreases, SiC will con- 
dense and coat the already existing carbon grains 
raising the relative emission of SiC. Finally, as 
mass-loss increases and the shell further expands 
and cools, MgS will condense last onto existing dust 
grains. In this scenario, one would expect no SiC 
emission for low M, then a gradually increasing SiC 
feature for intermediate mass-loss rates. Higher M 



would be characterized by a decrease in the SiC 
strength due to MgS coating the SiC dust grains. 

II. SiC -> C -> MgS: In this case, SiC would con- 

dense first at low mass-loss rates and is charac- 
terized by a sharp increase in SiC-to-C ratio. Car- 
bon condenses next as mass-loss progresses and the 
circumstellar envelope cools, reducing the SiC-to- 
C ratio. MgS forms last on the existing existing 
grains further reducing the SiC strength. 

III. SiC -> (C, SiC) -> MgS: This scenario is char- 

acterized by an initially rapid increase of SiC emis- 
sion. If C/O ratio, metallicity, and density are 
such that SiC and C have comparable condensation 
temperatures, then we would observe a subsequent 
SiC-to-C ratio that is constant with M until MgS 
begins to condense. Afterwards, the observed SiC 
strength would decline. 

lLagadec et al.l (|2007t ) argue for a SiC — > C conden- 
sation sequence for C-stars in the Galaxy. Figures fTTk 
and [TTb show that the relative emission of SiC-to-C for 
Galactic objects increases at low mass-loss rates; but 
due to the large scatter and high uncertainties at low 
M gas , we are unable to determine with absolute cer- 
tainty whether the trend levels off or continues to in- 
crease until 10 -6 M Q yr _1 . If it does indeed flatten 
at Mgas ~ 10~ 7 ' 3 M© yr _1 , this would suggest a sce- 
nario where SiC forms initially; then, as the mass-loss 
rate increases, C and SiC condense simultaneously. How- 
ever, based on the increasing trend in Figure [TTb as well 
as observational evidence of presolar grains indicating 
that G alactic carbon stars prod uce mainly uncoated SiC 
grains ([Bernatowicz et aDll987T ) , it is more likely that the 
SiC-to-C ratio increases until M gas ~ 10~ 6 M© yr _1 sup- 
porting aC-» SiC condensation sequence. The subse- 
quent decline in SiC strength correlates with an increase 
in MgS emission indicating that MgS is condensing onto 
the SiC grain. 
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quences at the top and plot the expected SiC-to-C ratio with mass- 
loss rate for each scenario. The dashed line indicates the onset of 
MgS grain condensation. 

Due to the lack of objects with low mass-loss rates ob- 
served in the Magellanic Clouds, it is difficult to discern 
whether SiC or graphite grains form initially. However, 
the fact that we observe a relatively large scatter in the 
SiC-to-C ratio for objects without MgS emission (espe- 
cially for the SMC sample) suggests near-simultaneous 
condensation of SiC and c arbon, supporting a Si C — * (C, 
SiC) sequence. Similarly, lLagadec et alj (|2007| ) suggest 
that SiC and C form close together in the LMC, but 
that SiC forms last in the SMC stars. A low initial Si 
abundance may also contribute to the observed low SiC 
strength in these environments, particularly, in the SMC. 
Finally, similar to the Galactic stars, a decline in SiC 
strength is observed as MgS increases for both the LMC 
and SMC objects. 

Based on the positions of the model from iWoitkd pflfih 
in the three panels of Figure [121 stars with higher metal- 
licities would more likely follow a C — ► SiC — > MgS 
condensation sequence in an outflowing wind. Lower- 
ing the photospheric temperature and/or the C/O ratios 
would shi ft the model atmosph ere curves down and to 
the right (|j0rgensen et al.| [l992). This would not affect 
the condensation sequence for Galactic carbon stars due 
to the initially large spread between the condensation 
zones for the different dust species; but it could poten- 
tially flip the C/SiC sequence shown for the Magellanic 
Clouds samples, especially for the SMC stars. However, 
the increased efficiency of carbon dredge-up leading to 
relatively higher C/O ratios at lower metallicities effec- 
tively reduces our ability to shift the LMC and SMC 
models. Thus, this supports the above suggestion that 
carbon stars in both Magellanic Clouds have SiC — ■> (C, 
SiC) — > MgS sequences. 

6. CONCLUSIONS 

We analyze a combined total of 107 carbon-rich AGB 
stars in the LMC, SMC, and Milky Way to investigate 
the effects of metallicity on dust and molecular species 



in these objects' circumstellar envelopes. Using modi- 
fied blackbody fits, we isolated and studied the spectral 
features throughout the 5-38 [iva. range including the P, 
Q, and R branches of C2H2 absorption at 7.5, 13.7, and 
14 /im; the 11.3 /im SiC emission feature; and the 30 /im 
MgS emission feature. 

The modified blackbody fits deviate from the Manch- 
ester method used in SZL by providing a more robust 
method of determining the continua above and below 
the aforementioned spectral features. In particular, this 
work improves on the existing IRS studies from SZL by 
allowing us to isolate the P and R branches in the 14 /im 
C2H2 band. For a given gas mass-loss rate, we find that 
LMC and SMC objects have comparable C2H2 equiva- 
lent widths, contrary to the findings of SZL. However, 
in agreement with SZL, we also find that carbon stars 
within the LMC and SMC have much larger C2H2 equiv- 
alent widths than Galactic carbon stars for a given M gas , 
which can be explained by a combination of more effi- 
cient carbon dredge-up and less dust dilution at lower 
metallicities. 

For the 13.7 fim C2H2 Q branch, we show that dust di- 
lution plays a prominent role on the shape and strength 
of this feature at lower mass-loss rates. However, at 
higher M, C2H2 absorption overcomes the effects of dust 
dilution due to increasing densities. Since veiling of the 
absorbing molecules by dust does not increase with the 
mass-loss rates, we c onfirm the previous results from 
iMatsuura et alj (|2006ft and SZL that the 13.7 fim C 2 H 2 
absorption band must originate throughout the circum- 
stellar dust shell. Conversely, the effects of dust veiling 
on the P and R branches become apparent at higher gas 



mass- loss rates {M gas > 10 



-6.5 



M 



yr" 



This sup- 



ports the idea that the optimal conditions for P and R 
branch transitions occur at depths below those favoring 
Q branch transitions. 

For objects with low to moderate mass-loss rates, we 
find that the molecular absorption band at 14 /jm af- 
fects the position of \sic- The adjoining P branch of 
strong acetylene absorption effectively reduces flux on 
the red side of the SiC feature causing Xsic to shift to- 
wards the blue. This effect occurs more prominently in 
the Magellanic Clouds because of the increased C2H2 
abundance due to enhanced carbon dredge-up. Simi- 
larly, an unidentified absorption feature at around 10 /im 
shifts Xsic towards longer wavelengths for objects with 
[6.4] — [9.3] < 0.2. SiC self-absorption dominates the 
shift in Asic for objects with large dust optical depths 
([6.4] — [9.3] color) and small C2H2 absorption at 14 /im. 
Differences in trends between Galactic, LMC, and SMC 
sources reflect the variation of C2H2 with metallicity for 
a given [6.4] — [9.3] color. 

Based on emission features of the SiC and MgS dust 
components with respect to the mass-loss rate, we have 
developed a method to determine the dust condensation 
history for a given metallicity system. We find evidence 
to support aC^ SiC condensation sequence in Galactic 
objects. For C-stars within the Magellanic Clouds, we 
suggest that a SiC — > (C, SiC) sequence best fits the ob- 
served data. Formation and increase of the MgS feature 
with Mg as corresponds to a decrease in the SiC strength 
indicating that MgS inhibits the SiC emission by coating 
the SiC grains. The condensation of MgS in circumstellar 



14 



Leisenring et al. 



atmosphere appears to be delayed until higher pressures 
for decreasing metallicities. 

The precision of much of this analysis depends on our 
ability to accurately determine individual gas mass-loss 
rates. In this work, we make three major assumptions 
that lead to uncertainties when deriving M gas : 1) the 
objects in each galaxy are comprised of the same metal- 
licity; 2) the dust-to-gas ratio scales with metallicity; 3) 
and M gas can be derived from the dust-to-gas ratio and 
Mdust- Any of these assumptions may contribute to the 
large scatter observed in Figure [Tib for the non-MgS 
bearing objects in the Magellanic Clouds as well as de- 
viations from the analysis by SZL when comparing gas 
mass-loss rates. Therefore, we are anticipating more ac- 
curate determinations of the gas mass-loss rate from di- 
rect measurements of molecular species such as CO, par- 
ticularly in the Magellanic Clouds. Differences between 
CO-determined M gas and those derived in this and pre- 
vious works will further improve our understanding of 
dust formation efficiencies as a function of metallicity. 
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07 








MSX LMC 1209 





87±0.01 





259±0 


006 





232±0 


011 





074±0 


007 





124±0 


004 


11 


12±0 


01 








MSX LMC 1213 





77±0.01 





127±0 


003 





114±0 


015 





043±0 


013 





159±0 


003 


11 


09±0 


03 





481±0 


011 


MSX LMC 1298 


o 


61±0.01 





066±0 


002 








o 


050±0 


003 


o 


106±0 


001 


11 


39±0 


02 


o 


083±0 


003 


MSX LMC 1308 





57±0.01 





181±0 


006 





609±0 


013 





041±0 


008 





037±0 


007 


10 


07±0 


11 








MSX LMC 1383 


1 


26±0.01 





176±0 


004 





262±0 


007 





026±0 


004 





052±0 


002 


11 


21±0 


07 





633±0 


005 


TV TCI "V-" TUT /~1 r* O A 

MSX LMC 634 





60±0.01 





078±0 


004 





118±0 


012 





055±0 


007 





044±0 


003 


11 


34±0 


08 








MSX LMC 768 





85±0.01 





122±0 


002 





191±0 


005 





042±0 


003 





085±0 


002 


11 


27±0 


03 





072±0 


005 


MSX LMC 774 


1 


08±0.01 





226±0 


004 





273±0 


005 





056±0 


003 





075±0 


002 


11 


11±0 


06 





349±0 


006 


MSX LMC 783 





68±0.01 





114±0 


003 





271±0 


010 





048±0 


008 





128±0 


003 


11 


13±0 


04 








MSX LMC 787 





96±0.01 





195±0 


004 





327±0 


011 





023±0 


006 





110±0 


003 


11 


10±0 


01 





645±0 


096 


MSX LMC 841 





68±0.01 





141±0 


004 





094±0 


008 





047±0 


006 





081±0 


003 


11 


06±0 


06 








MSX LMC 937 





91±0.01 





130±0 


002 





109±0 


005 





058±0 


003 





101±0 


002 


11 


23±0 


03 





278±0 


004 


MSX LMC 950 





88±0.01 





147±0 


004 





361±0 


009 





075±0 


006 





112±0 


004 


11 


10±0 


05 





287±0 


017 


MSX LMC 971 


1 


12±0.01 





166±0 


005 





433±0 


007 





040±0 


005 





069±0 


002 


11 


19±0 


01 





422±0 


006 


MSX LMC 974 





74±0.01 





055±0 


005 





223±0 


009 





063±0 


006 





049±0 


003 


11 


01±0 


08 








Ziilstra et al. (2006) 


Sample: 






































GRRV 38 





24±0.03 





384±0 


017 





237±0 


032 





135±0 


023 





173±0 


016 


11 


17±0 


08 








IRAS 04557-6753 


1 


07±0.01 





218±0 


006 





195±0 


007 





060±0 


003 





091±0 


004 


11 


21±0 


07 





461±0 


004 


IRAS 05009-6616 





90±0.00 





084±0 


002 





111±0 


004 





051±0 


002 





122±0 


001 


11 


21±0 


02 





267±0 


003 


IRAS 05112-6755 


1 


12±0.01 





180±0 


003 





386±0 


004 





058±0 


002 





052±0 


001 


11 


08±0 


01 





317±0 


002 


IRAS 05113-6739 


1 


05±0.02 





187±0 


004 





085±0 


007 





044±0 


005 





066±0 


002 


11 


23±0 


08 





299±0 


004 


IRAS 05132-6941 


1 


05±0.01 





154±0 


004 





219±0 


010 





062±0 


007 





094±0 


006 


11 


24±0 


08 





138±0 


006 


IRAS 05190-6748 


1 


25±0.01 





172±0 


003 





251±0 


006 





063±0 


004 





044±0 


005 


11 


15±0 


20 





418±0 


001 


IRAS 05278-6942 


1 


11±0.01 





135±0 


004 





295±0 


009 





074±0 


007 





067±0 


005 


11 


02±0 


12 





502±0 


002 


IRAS 05295-7121 





95±0.01 





108±0 


004 





139±0 


007 





063±0 


004 





066±0 


003 


11 


20±0 


07 





169±0 


003 


IRAS 05360-6648 


1 


06±0.01 





157±0 


004 





279±0 


007 





042±0 


004 





125±0 


003 


11 


07±0 


01 





691±0 


005 


MSX LMC 219 


1 


13±0.01 





256±0 


003 





541±0 


006 





048±0 


004 





048±0 


002 


11 


18±0 


05 





384±0 


004 


MSX LMC 341 





96±0.01 





142±0 


006 





166±0 


007 





050±0 


003 





124±0 


002 


11 


01±0 


03 





672±0 


008 


MSX LMC 349 





qo_i_n ni 





113±0 


003 





322±0 


007 


o 


045±0 


003 


o 


122±0 


003 


11 


04±0 


01 


o 


740±0 


008 


MSX LMC 441 


1 


12±0.01 





195±0 


005 





407±0 


007 





046±0 


004 





082±0 


002 


11 


07±0 


03 





647±0 


004 


MSX LMC 443 





89±0.01 





261±0 


004 





364±0 


009 





037±0 


006 





119±0 


003 


11 


11±0 


03 





534±0 


004 


MSX LMC 494 





55±0.01 





186±0 


003 





065±0 


004 





036±0 


003 





095±0 


002 


11 


23±0 


03 








MSX LMC 601 





89±0.01 





035±0 


003 





109±0 


005 





055±0 


003 





043±0 


003 


11 


11±0 


11 








MSX LMC 679 





89±0.01 





161±0 


003 





267±0 


008 





035±0 


005 





131±0 


005 


11 


07±0 


07 





534±0 


003 


MSX LMC 743 


1 


03±0.01 





279±0 


004 





425±0 


007 





028±0 


004 





107±0 


003 


11 


11±0 


04 





701±0 


003 


MSX LMC 749 





60±0.01 





146±0 


003 





320±0 


008 





067±0 


005 





063±0 


006 


11 


02±0 


15 








MSX LMC 754 





60±0.02 





122±0 


004 





554±0 


012 





013±0 


008 





132±0 


004 


11 


10±0 


01 








MSX LMC 967 





59±0.01 





055±0 


004 





168±0 


012 





062±0 


011 





065±0 


006 


11 


23±0 


21 








NGC 1978 IR1 





55±0.01 





122±0 


002 





383±0 


007 





057±0 


003 





088±0 


002 


11 


08±0 


05 








TRM 72 





78±0.01 





051±0 


004 





121±0 


007 





069±0 


004 





064±0 


002 


11 


18±0 


05 








TRM 88 





76±0.02 





121±0 


006 





178±0 


006 





055±0 


003 





058±0 


007 


10 


97±0 


13 








Sloan et al. (2006) Sample: 

MSX SMC 033 0.51±0.01 





052±0 


002 





081±0 


006 





050±0 


003 





059±0 


002 


10 


96±0 


06 








MSX SMC 036 





76±0.01 





299±0 


005 





401±0 


008 





054±0 


005 





042±0 


002 


11 


06±0 


04 








MSX SMC 044 





53±0.01 





041±0 


003 











066±0 


008 





079±0 


003 


10 


98±0 


07 








MSX SMC 054 





75±0.01 





236±0 


002 





209±0 


006 





057±0 


004 





139±0 


001 


11 


19±0 


01 





177±0 


020 


MSX SMC 060 





96±0.00 





215±0 


003 











037±0 


003 





030±0 


002 


11 


27±0 


08 





228±0 


027 


MSX SMC 062 





65±0.02 





129±0 


006 





055±0 


012 





053±0 


007 





055±0 


003 


11 


13±0 


09 








MSX SMC 066 





43±0.01 





091±0 


002 





464±0 


012 


o 


040±0 


004 


o 


025±0 


002 


11 


24±0 


01 








MSX SMC 091 





65±0.02 





358±0 


010 





329±0 


022 





059±0 


017 





121±0 


007 


11 


23±0 


07 








MSX SMC 093 





46±0.01 





241±0 


004 


1 


077±0 


014 





027±0 


008 





016±0 


003 


10 


81±0 


11 








MSX SMC 105 





87±0.01 





297±0 


004 





085±0 


006 





029±0 


004 





120±0 


003 


11 


12±0 


1 





487±0 


014 


MSX SMC 142 





29±0.03 





165±0 


007 











090±0 


030 





131±0 


016 


11 


56±0 


28 








MSX SMC 159 





87±0.01 





204±0 


002 





285±0 


005 





047±0 


004 





124±0 


001 


11 


04±0 


02 





417±0 


010 


MSX SMC 162 





50±0.01 





091±0 


004 





393±0 


016 





006±0 


011 





078±0 


004 


11 


09±0 


10 








MSX SMC 163 





68±0.01 





171±0 


003 





276±0 


006 





054±0 


003 





119±0 


002 


11 


15±0 


02 





236±0 


008 


MSX SMC 198 





62±0.01 





141±0 


003 





009±0 


011 





073±0 


005 





055±0 


002 


11 


10±0 


05 








MSX SMC 200 





42±0.01 





007±0 


003 





095±0 


007 





069±0 


005 




















MSX SMC 202 





45±0.02 





158±0 


006 





078±0 


019 





002±0 


017 





037±0 


007 


10 


65±0 


23 








MSX SMC 209 





70±0.01 





198±0 


002 





350±0 


007 





047±0 


003 




















MSX SMC 232 





63±0.01 





100±0 


002 





270±0 


009 





061±0 


004 





009±0 


002 


10 


92±0 


15 








Lagadec et al. ( 

GM780 


2007) Sample: 

0.41±0.01 





340±0 


010 





164±0 


024 





035±0 


016 





246±0 


007 


11 


07±0 


05 








IRAS 00554-7351 


0.94±0.01 





145±0 


005 





186±0 


009 





078±0 


007 





085±0 


005 


11 


27±0 


09 





291±0 


004 


ISO 00548 


0.41±0.01 





060±0 


003 





023±0 


008 





064±0 


005 





086±0 


002 


11 


30±0 


03 








ISO 00573 


0.30±0.01 





167±0 


006 





197±0 


013 





078±0 


008 





007±0 


007 


11 


72±0 


56 








ISO 01019 


0.41±0.02 





137±0 


004 





001±0 


012 





039±0 


009 





056±0 


005 


11 


25±0 


10 









Metallicity Effects on Carbon Stars 

TABLE 2 — Continued 
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[6.4] - [9.3] 




EW(7.5) 




EW(14) 


EW(13.7) 








XsiC 










Name 


(mag) 




(/^m) 




(/mi) 




H 




L/C (SiC) 




(/im) 




L/C (MgS) 


LEGC 105 


0.36±0.02 





094±0.002 


0. 


198zb0. 


010 


n n c a i n 
U.Uo4±U. 


005 


0. 


044±0.003 


11 


09±0. 


07 








NGC 419 IR1 


0.49±0.01 





152±0.002 


0. 


766zb0. 


.008 


U.U41±U. 


.004 


















NGC 419 LE 16 


0.34±0.02 





090±0.004 


0. 


185zb0. 


016 


n n^^J_n 
U.Ubb±U. 


012 





071±0.004 


11 


23±0. 


06 








NGC 419 MIR1 


1.36±0.01 





,278±0.005 


0. 


oOOdzt). 


007 


n n a o _i_n 


004 


0. 


023±0.003 


11 


,40±0. 


13 


0. 


409±0. 


003 


RAW 960 


0.23±0.01 





226±0.004 


0. 


1/11 i n 




n r\p, \ i n 
U.Uu4±U. 


nnn 

uuy 


0. 


063±0.005 


11 


,17±0. 


09 








Sloan et al. (2003) Sample: 
































AFGL 2392 


0.58±0.00 





,081±0.001 


0. 


028zb0. 


002 


n noo I n 
U.UzZ±U. 


001 


0. 


,312±0.000 


11 


.19±0. 


00 


0. 


336±0. 


002 


AFGL 3099 


1.21±0.00 





063±0.001 


0. 


008i0. 


.001 


n r\a. *7 i n 


.001 





130±0.000 


11 


,43±0. 


00 





,427±0. 


001 


AFGL 341 


1.67±0.01 





080±0.001 


0. 


c\a o_i_n 
U4ZztU. 


UU1 


n ncn I n 


nm 
UU1 


0. 


032±0.000 


12 


,12±0. 


01 


0. 


318±0. 


001 


AFGL 940 


0.77±0.00 





040±0.001 


0. 


OlOzbO. 


003 


n n a 1 i n 
U.U41±U. 


002 


0. 


176±0.001 


11 


,27±0. 


01 


0. 


311±0. 


.002 


C 2178 


0.42±0.00 





155±0.001 


0. 




f\C\ I 

UU4 


n n/i 1 i n 
U.U4I±U 


nno 
UUz 





200±0.001 


11 


23±0. 


01 








C 2429 


0.19±0.00 





166±0.001 








n nQO_Ln 


nno 
UUz 


0. 


614±0.001 


11 


36±0. 


.00 








Case 181 


0.61±0.00 





047±0.001 


0. 


081±0 


.001 


n n c o _Ln 


.001 


0. 


199±0.000 


11 


,25±0. 


.00 


0. 


186±0. 


.000 


DO 16793 


1.29±0.00 





043±0.000 


0. 


252-J-.0 


.001 


n n a o_i_n 


.001 





078±0.000 


11 


30±0. 


01 





,454±0. 


.000 


DO 40123 


0.29±0.00 





311±0.001 


0. 


lOliO 


.002 


n nt;7-Ln 


.001 





275±0.000 


11 


,08±0. 


00 








HV Cas 


0.35±0.00 





146±0.000 


0. 


c\c\a i n 
UU4±U 


UUz 


n n/i q i n 
U.U4oztU 


nm 
UU1 


0. 


216±0.000 


11 


.34±0. 


.00 








IRAS 21489+5301 


1.36±0.01 





102±0.001 


0. 


072±0 


.002 


n n /i /i _Ln 
U.U44±U 


.002 


0. 


018±0.001 


12 


09±0. 


02 


0. 


347±0. 


.001 


IRC+00365 


0.67±0.00 





077±0.000 




' ' ' 




U.Uo2±U 


.001 


0. 


178±0.000 


11 


32±0. 


00 





206±0. 


.000 


IRC+20326 


0.89±0.00 





009±0.001 


0. 


030±0 


.001 


n no '"> I n 
U.UUo±U 


.000 


0. 


078±0.000 


11 


.44±0. 


01 


0. 


354±0. 


.000 


IRC+40540 


1.16±0.00 





064±0.000 


0. 


OOliO 


.000 


n no/? I n 


.000 


0. 


148±0.000 


11 


33±0. 


.00 


0. 


310±0. 


.000 


IRC+50096 


0.68±0.00 





111±0.001 


0. 


Z4uztU 


UU1 


n n*3 1 i n 


nm 
UU1 





257±0.000 


11 


13±0. 


00 





,144±0. 


.000 


IRC-10095 


0.16±0.00 





292±0.001 




' ' ' 




n n 1 q I n 
U.Ulo±U 


.002 





140±0.001 


11 


35±0. 


01 








IRC-10122 


0.38±0.00 


0. 


093±0.000 


0. 


140zb0 


.003 


n c\q i i n 
U.Uol±U 


.001 


0. 


243±0.000 


11 


,15±0. 


00 








RU Vir 


0.51±0.00 





076±0.001 


0. 


no c i n 


C\f\ "1 

UU1 


n m 1 i n 


nm 
UU1 


0. 


320±0.000 


11 


,08±0. 


00 








RY Dra 


O.lliO.OO 





115±0.000 




' ! ' 




n nc 1 i n 


nm 
UU1 


0. 


099±0.000 


11 


61±0. 


00 








R For 


0.36±0.00 





142±0.001 


0. 


073zb0 


.002 


n no *7 I n 
U.Uzo±U 


.001 


0. 


264±0.000 


11 


,13±0. 


.00 








R Scl 


0.27±0.00 





390±0.001 


0. 


712±0 


.003 


n A77±n 
U.U/ r ±U 


.001 


0. 


154±0.001 


11 


09±0. 


.00 








SS Vir 


0.36±0.00 





420±0.000 


0. 


274i0 


.009 


n n cc o i n 


.004 





317±0.001 


11 


,17±0. 


.00 








S Cep 


0.30±0.00 





160±0.000 








n n/i c i n 


n n i 
UU1 


0. 


253±0.001 


11 


,15±0. 


01 








S Set 


0.14±0.00 





293±0.001 


0. 


UzZztU 




n nQC-un 


nno 


0. 


160±0.000 


11 


,15±0. 


01 








TU Tau 


0.08±0.00 












0.128±0 


.007 


0. 


270±0.001 


11 


,44±0. 


01 








T Dra 


0.51±0.00 





162±0.001 


0. 


137±0 


.003 


0.026±0 


.001 





344±0.000 


11 


,07±0. 


.00 


0. 


149±0. 


002 


U Cam 


O.lliO.OO 





065±0.000 


0. 


095±0 


.002 


0.075±0 


.001 


0. 


279±0.000 


11 


.34±0. 


.00 








V460 Cyg 


0.16±0.00 





125±0.000 


0. 


067=b0 


.004 


0.077±0 


.003 





049±0.001 


11 


56±0. 


02 








VX And 


0.12±0.00 





214±0.000 


0. 


260±0 


.004 


O.lOliO 


.002 





115±0.000 


11 


39±0. 


.00 








V Aql 


O.lliO.OO 





146±0.001 


0. 


156±0 


.002 


0.097±0 


.001 


0. 


098±0.000 


11 


38±0. 


.00 








V CrB 


0.40±0.00 





166±0.001 


0. 


180±0 


.003 


0.053±0 


.001 


0. 


387±0.000 


11 


,12±0. 


.00 


0. 


140±0. 


.005 


VCyg 


0.52±0.00 





196±0.000 


0. 


351=b0 


.003 


0.036±0 


.001 





309±0.000 


11 


01±0. 


.00 


0. 


060±0. 


.001 


W Ori 


0.06±0.00 





058±0.000 


0. 


006±0 


.001 


0.082±0 


.001 


0. 


144±0.000 


11 


,47±0. 


.00 








Y CVn 


0.15±0.00 





179±0.001 


0. 


029±0 


.002 


0.067±0 


.002 


0. 


068±0.000 


11 


50±0. 


.01 









Note. — For all objects in each sample, we show the [6.4] — [9.3] color as well as strengths of the molecular and dust features. The molecular 
features arc measured in terms of equivalent widths (EW) in microns while the dust emission bands show a linc-to-continuum ratio (L/C). We 
also give the central wavelength of the SiC feature (Xsic)- 
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TABLE 3 
Mass-loss Rates 



Name 


logM,;,, ,, a lo 
fMn vr - 


g Mgas b 

i) 

) 


Jylrlii a-h 

tt Lfc o L 

Ref. 


Type c 
(Ref.) 


Name 


log Af^ T , ^ a log 
fMn vr -1 


Mgas b 


Adrln, q + 
(j, LL j L 

Ref. 


Type c 
(Ref.) 


Our Sample 




















TR A Cl C\AAQ« fiQ^S 

iri/\o u^^yo-oyoo 


7 ^ 


-o.zo 


i 




M^"V T Mr 1 7fic 
IVloyv i^lvlvy I Oo 


7 ^ 


^ 9^ 


i 




LI-LMC 603 


-6.76 


-4.46 


i 




MSX LMC 774 


-7.25 


-4.95 


i 




LI-LMC 1817 


-7.06 


-4.76 


i 




MSX LMC 783 


-7.80 


-5.50 


i 




MSX LMC 1205 


-7.81 


-5.51 


i 




MSX LMC 787 


-7.40 


-5.10 


i 




MSX LMC 1209 


-7.52 


-5.22 


i 




MSX LMC 841 


-7.80 


-5.50 


i 




MSX LMC 1213 


-7.67 


-5.37 


i 




MSX LMC 937 


-7.47 


-5.17 


i 




MSX LMC 1298 


-7.91 


-5.61 


i 




MSX LMC 950 


-7.51 


-5.21 


i 




MSX LMC 1308 


-7.98 


-5.68 


i 




MSX LMC 971 


-7.20 


-4.90 


i 




MSX LMC 1383 


-7.04 


-4.74 


i 




MSX LMC 974 


-7.71 


-5.41 


i 




MSX LMC 634 


-7.93 


-5.62 


i 














Ziilstra et al. (2006) Sample: 


















GRRV 38 


-8.64 


-6.34 


i 


R(2) 


MSX LMC 441 


-7.26 


-4.96 


2 


R(2) 


IRAS 04557-6753 


-7.22 


-4.92 


2 


N (2) 


MSX LMC 443 


-7.60 


-5.30 


2 


R(2) 


IRAS 05009-6616 


-7.33 


-5.03 


2 


N (2) 


MSX LMC 494 


-8.15 


-5.85 


2 


R 2 


IRAS 05112-6755 


-7.05 


-4 74 


2 


IN \ J.) 


T A/TP fifll 


7 

- 1 .oy 


^ 9Q 


2 


it (Z) 


IRAS 05113-6739 


-7.22 


-4.92 


2 


N (2) 


MSX LMC 679 


-7.54 


-5.24 


2 


R(2) 


IRAS 05132-6941 


-7.32 


-5.02 


2 


R(2) 


MSX LMC 743 


-7.43 


-5.13 


2 


R(2) 


IRAS 05190-6748 


-6.94 


-4.64 


2 


N (2) 


MSX LMC 749 


-7.70 


-5.40 


2 


R(2) 


IRAS 05278-6942 


-6.90 


-4.60 


2 


R(2) 


MSX LMC 754 


-8.07 


-5.77 


2 


R(2) 


IRAS 05295-7121 


-7.35 


-5.05 


2 


N (2) 


MSX LMC 967 


-7.80 


-5.49 


2 


R(2) 


IRAS 05360-6648 


-7.39 


-5.09 


2 


N (2) 


NGC 1978 IR1 


-7.94 


-5.64 


2 


N (2) 


MSX LMC 219 


-7.30 


-5.00 


2 


R(2) 


TRM 72 


-7.54 


-5.24 


2 


N (2) 


MSX LMC 341 


-7.43 


-5.13 


2 


R(2) 


TRM 88 


-7.78 


-5.48 


2 


R(2) 


MSX LMC 349 


-7.40 


-5.10 


2 


R(2) 










Sloan et al. (2006 


Sample: 


















MSX SMC 033 


-7.65 


A 


9 


R (0\ 
it (Z) 


1\/TQV QA/Tr 1 1/10 


Q x.7 
-o.O / 


c; Q7 


Z 




MSX SMC 036 


-7.76 


-5.06 


2 


N (2) 


MSX SMC 159 


-7.62 


-4.92 


2 


N (2) 


MSX SMC 044 


-7.70 


-5.00 


2 


R(2) 


MSX SMC 162 


-8.19 


-5.49 


2 


N (2) 


MSX SMC 054 


-7.74 


-5.05 


2 


MSX SMC 163 


-7.78 


-5.08 


2 


R(2) 


MSX SMC 060 


-7.55 


-4.85 


2 


N (2) 


MSX SMC 198 


-8.12 


-5.43 


2 


R(2) 


MSX SMC 062 


-7.80 


-5.10 


2 


R(2) 


MSX SMC 200 


-7.89 


-5.19 


2 


MSX SMC 066 


-7.77 


-5.07 


2 


N (2) 


MSX SMC 202 


-8.36 


-5.66 


2 




MSX SMC 091 


-8.12 




2 


M^"V ^MP 9flQ 


7 7 A 


-o.uo 


n 


R (9\ 
it (Z) 


MSX SMC 093 


-8.60 


-5.90 


2 




MSX SMC 232 


-7.84 


-5.14 


2 


N (2) 


MSX SMC 105 


-7.67 


-4.97 


2 


r"(2) 










Lagadec et al. (2007) Sample: 


















GM780 


-8.22 


-5.52 


2 


R(2) 


LEGC 105 


-8.26 


-5.56 


2 


R (2) 


IRAS 00554-7351 


-7.19 


-4.49 


2 


N (2) 


NGC 419 IR1 


-8.00 


-5.30 


2 


N (2) 


ISO 00548 


-7.96 


-5.26 


2 


R(2) 


NGC 419 LE 16 


-8.59 


-5.89 


2 


N (2) 


ISO 00573 


-8.51 


-5.81 


2 


R(2) 


NGC 419 MIR1 


-7.07 


-4.37 


2 


R (2) 


ISO 01019 


-8.26 


-5.56 


2 


R(2) 


RAW 960 


-8.76 


-6.06 


2 


R (2) 


Sloan et al. (2003 


Sample: 


















AFGL 2392 


-7.74 


-5.74 


3 


M (5) 


RU Vir 


-7.96 


-5.96 


4 


M (7) 


AFGL 3099 


-7.10 


-5.10 


3 


M (6) 


RY Dra 


-9.82 


-7.82 


4 


S (7) 


AFGL 341 


-7.02 


-5.02 


4 


M (6) 


R For 


-8.10 


-6.10 


3 


M (7) 


AFGL 940 


-7.49 


-5.49 


4 


R Scl 


-9.15 


-7.15 


3 


S (7) 


C 2178 


-7.95 


-5.95 


4 




SS Vir 


-8.74 


-6.74 


4 


S (7) 


C 2429 


-9.11 


-7.11 


4 


I '(7) 


S Cep 


-8.71 


-6.71 


4 


M (7) 


Case 181 


-7.32 


-5.32 


4 


I (6) 


S Set 


-10.08 


-8.08 


4 


S (7) 


DO 16793 


-12.06 


-10.06 


4 


TU Tau 


-9.85 


-7.85 


4 


S (7) 


DO 40123 


-8.58 


-6.58 


4 


M (6) 


T Dra 


-8.41 


-6.41 


4 


M (7) 


HV Cas 


-8.05 


-6.05 


4 


M (5) 


U Cam 


-9.75 


-7.75 


4 


S (7) 


IRAS 21489+5301 


-6.80 


-4.80 


4 


V460 Cyg 


-10.51 


-8.51 


4 


S (7) 


IRC+00365 


-7.50 


-5.50 


4 


M (6) 


VX And 


-9.33 


-7.33 


4 


S (7) 


IRC+20326 


-7.27 


-5.27 


4 




V Aql 


-9.35 


-7.35 


4 


S (7) 


IRC+40540 


-7.08 


-5.08 


4 


M (7) 


V CrB 


-9.08 


-7.08 


4 


M (7) 


IRC+50096 


-7.81 


-5.81 


4 


M (7) 


VCyg 


-8.16 


-6.16 


4 


M (7) 


IRC- 10095 


-8.52 


-6.52 


4 




W Ori 


-9.87 


-7.87 


4 


S (7) 


IRC-10122 


-7.92 


-5.92 


4 


M (6) 


Y CVn 


-9.92 


-7.92 


4 


S (7) 



References. — (1) Relatio n fr om Ivan Loonl d2007l) : (2 ) IGroenewegen et "all (120071) : (3) iLe BerTrel U997T : (4) Relation from lLe Bertrd (119971 ): 
CSl lGroeixewegen et all d2002D : (6) \Louv> endr(ll993r i: r7l fBergeat &: Chevallieri TftOOSi ) 

a b 
Dust mass-loss rates are cither results of dust radiative transfer models of a given source or taken from relations based off" of these models. Gas 

mass-loss rates were derived using Equations [3] through [5] C Objects in the LMC and SMC exhibit Mira-likc pulsations: R stand for regular 

pulsation and N stands for non-regular pulsation. For Galactic objects, M stand for Mira, I stand for Irregular, and S stand for Seiniregular. 



